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Actuating periodically an elastic filament in a viscous liquid generally breaks the constraints of 
Purcell's scallop theorem, resulting in the generation of a net propulsive force. This observation 

■ suggests a method to design simple swimming devices - which we call "elastic swimmers" - where the 
actuation mechanism is embedded in a solid body and the resulting swimmer is free to move. In this 
paper, we study theoretically the kinematics of elastic swimming. After discussing the basic physical 
picture of the phenomenon and the expected scaling relationships, we derive analytically the elastic 
swimming velocities in the limit of small actuation amplitude. The emphasis is on the coupling 
between the two unknowns of the problems - namely the shape of the elastic filament and the 
swimming kinematics - which have to be solved simultaneously. We then compute the performance 
of the resulting swimming device, and its dependance on geometry. The optimal actuation frequency 
and body shapes are derived and a discussion of filament shapes and internal torques is presented. 
Swimming using multiple elastic filaments is discussed, and simple strategies are presented which 
result in straight swimming trajectories. Finally, we compare the performance of elastic swimming 

> . with that of swimming microorganisms. 

■ 

S: 

■ ■ I. INTRODUCTION 

■ 

C ■ 

The fluid mechanics of microorganism locomotion, pioneered more than fifty years ago by G.I. Taylor, has become 
, ^i - one of the most successful branches of biomechanics, with success in both the basic physical understanding of flow 
behavior and the quantitative prediction of kinematics and energetics of locomotion [1-8] . 

Recent technical advances have led to ever more precise fabrication at small scales (microns or less), prompting both 
theorists [9-13] and experimentalists [14] to design and analyze a series of simple low-Reynolds number swimmers. 
The experiment of Dreyfus et al. [14], in particular, reported locomotion in a sperm-like micro- swimmer, composed of 
a cargo (red blood cell) and a slender flexible filament made of a series of paramagnetic beads. In that case, actuation 
by oscillating transverse magnetic fields led to the generation of bending waves propagating along the filament and 
resulted in the motion of the micro- swimmer. In this system, the right-left symmetry was broken by the presence of a 
cargo and led to a preferential tip-to-base propagation of the bending waves, resulting in locomotion in the direction 

■ base-to-tip. 

An alternative way to break the symmetry in a similar system would be to build-in the asymmetry in the actuation. 
In particular, if an elastic filament is periodically actuated at one of its extremities in a viscous liquid, the resulting 
motion will lead to the propagation of bending waves and, in general, propulsive forces. This idea was originally 
proposed by Edward Purcell [8]. Physically, actuating an elastic filament allows one to break the constraints of the 
"scallop theorem" - which states that a body performing a reciprocal motion at low Reynolds number cannot propel 
itself - by allowing the boundary conditions on the fluid problem - that is, the shape of the filament - to be itself a 
function of the fluid flow. The original theoretical study on this problem was proposed by Wiggins and Goldstein [15], 
who showed that the amplitude of the actuated elastic filament satisfies a hyperdiffusion equation. This equation was 
also derived in earlier work by Machin in the context of wave propagation in the flagella of swimming microorganisms 
[16, 17]. A similar theoretical treatment was proposed as a simple model of the sliding filament model of eukaryotic 
axonemal beating by coupling the elasto-hydrodynamics problem with models for the behavior of active molecular 
motors [18, 19]. 

The main features of this problem have been successfully exploited experimentally to measure the bending modulus 
of biopolymers (actin filaments and microtubules), either using thermal fluctuations [20] or using an active actuation 
[21, 22]. Related studies include the dynamics of magnetic filaments [23-25], the three-dimensional actuation and 
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instabilities of flexible filaments [26-29] and the exploitation of symmetry-breaking to pump fluid in a channel [30]. 

In this paper, we consider the case where the actuated flexible filaments are exploited for locomotion purposes. We 
consider a prototypical micro- swimmer composed of a solid body and an elastic slender filament (see Fig. 1). The 
filament is fixed to the body and its base- angle is varied sinusoidally at frequency uj, in two or three dimensions, by 
a mechanism embedded in the swimmer body. This generates the propagation of bending waves down the elastic 
filament and propels the swimmer forward. This design for an "elastic swimmer" is the simplest locally-forced low- 
Reynolds swimmer exploiting the interplay of fluid drag and bending rigidity for propulsive purposes, and as a result, 
its swimming kinematics and energetics are of fundamental interest. 

As we will see below, the locomotion of the elastic swimmer also turns out to be an interesting mathematical 
problem. Indeed, to characterize the swimmer completely, two main problems need to be solved for, namely (1) 
the periodic shape of the elastic filament and (2) the kinematics of swimming. However, these two problems cannot 
be solved independently. Along the filament, drag forces and bending forces balance. On one hand, the swimming 
kinematics affects the drag and therefore the filament shape. On the other hand, the shape influences the viscous 
propulsive force and therefore the overall swimming velocity. As a result, these two problems have to be solved 
simultaneously, a fact which results - as will be seen below - in the appearance of integro-differential equations. This 
feature might have been overlooked by previous analytical studies. 

Numerical simulations of the elastic swimming problem were presented by Lagomarsino et al. [31] using particle- 
based methods (see also Ref. [32]). However, in this study, the filament was actuated by external forces and torques, 
and as a result is fundamentally different from the self-contained force-free and torque-free swimmer which we consider 
in this paper. Moreover, in the case of small amplitude oscillations of the actuation point, the simulations by 
Lagomarsino et al. obtain a constant swimming velocity for long filaments (long in the sense L ^> i^, see below), 
whereas in fact the velocity should decrease to zero because of excessive drag on the filament. This discrepancy is 
resolved in our paper. 

Numerical simulations of the three-dimensional actuation (rotation) of the filament were presented by Manghi et al. 
[26] using particle-based methods which include hydrodynamic interactions (similar to those used to study polymer 
dynamics). However, the simulations by Manghi et al. obtain swimming even in the case where the body sizes shrink 
to zero, a fact which also violates torque balance for a torque- free swimmer at zero Reynolds number [33, 34]. 

Recently, Yu et al. has performed a macro-scale experiment aimed at measuring the propulsive force generated by 
actuated elastic filaments in Stokes flows and compared it with existing theories [35]. The filaments were fixed in 
space with base-angles which were actuated sinusoidally and very good quantitative agreement was found between 
the measured propulsive force and that predicted by the small-amplitude theoretical study of Wiggins and Goldstein 
[15]. 

This paper is organized as follows. In §11 we present the basic physical picture for the coupling of hydrodynamics 
and bending forces in actuated filaments. We estimate the optimal actuation conditions of the filament and derive 
the expected scalings for the swimming speed of the swimmers. In §111 we derive the swimming kinematics of the 
elastic swimmers analytically in the limit of small actuation amplitude. The assumptions necessary to perform the 
calculation are clearly stated, and the final results are six analytical formulae for the three-dimensional trajectory of 
the swimmer (Eq. 49). The performance of the elastic swimmer is discussed in §IV. In particular, we characterize 
optimal swimmers as well as the difference between filament shapes for free-swimming versus fixed actuated filaments. 
Elastic swimming with more than one flexible filament is discussed in §V and we show that steady swimming on a 
straight line can be obtained with six filaments. Finally, a discussion of the results and a comparison of the swimmer 
performance with swimming microorganisms are presented in §VI. 

II. PHYSICAL PICTURE 
A. Elasto- hydrodynamics 

As explained above, we consider in this paper the prototypical elastic swimmer displayed in Fig. 1. We denote by 
L the length of the filament, r its radius, A its bending stiffness, and ^± its normal drag coefficient i.e. the viscous 
force exerted by the fluid per unit length of the filament for motion perpendicular to its length [3-5]. Let us first 
consider the case where the filament is actuated but not free to move [15, 35] . If y denotes the typical amplitude of 
a material point at a distance x along the filament, the balance between local viscous drag and bending forces on the 
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FIG. 1: Notations for the elastic swimmer. The filament has length L and radius r. The distance between the center of mass 
of the body and the actuation point is denoted a. The frame {e x ,e y ,e z } moves with the swimmer body whereas the frame 
{ei, e2, 63} is fixed. 



filament results in a hyperdiffusion equation for small- amplitude motion, 

dy _ d A y 

For a given actuation frequency lj, inspection of Eq. (1) shows the appearance of an intrinsic bending- hydrodynamics 
length scale, £ u = (A/uj^) 1 ^ [15]. If L <C bending forces win and the filament is straight. On the contrary, if 
L^> £ L0 drag forces win and the portions of the filament located at a distance larger than i u from the actuation point 
are essentially straight, a feature which is expressed mathematically by an exponentially decay of the amplitude of 
the solution to the elasto-hydrodynamics problem over the length scale t u [15]. 

In the case where the actuated filament is embedded in a swimming device, the limit L <C £ u results in a reciprocal 
actuation of the filament, and therefore - by use of the scallop theorem - the swimming velocity is zero [8]. The 
other limit L^> £^ leads to a constant value of the propulsive force, as there are no contribution from any portion of 
the filament beyond L ~ f^, and since in this case the viscous drag is large, and this limit also results in swimming 
velocity decreasing to zero. As a consequence, we expect that the optimal swimming will be obtained for a filament 
length L ~ £ u , taking full advantage of the drag-induced bending of the filament while keeping the overall drag on 
the swimmer low [31, 32]. This result is confirmed in §IV where we compute the optimal elastic swimmers. 



B. Optimal forcing 

Before considering the swimming kinematics, we show here using scaling arguments that it is energetically favorable, 
when actuating the elastic filament periodically, to only use the fundamental frequency. Let us call T the actuation 
period and e its amplitude. Such periodic actuation is necessary in order to constantly generate bending deformation 
and therefore propulsion. The mechanical forcing at the base of the filament can, in principle, include the fundamental 
frequency 27r/T, and all other harmonics. For small-amplitude motion, the propulsive force is quadratic in the filament 
dynamics, so by orthogonality we can study each frequency independently. Let us consider forcing at a given frequency 
uo. The elastic propulsive force scales as [15, 35] 



L ^d X (2) 

dxdx* dx - (2) 



For an optimal propulsion, the filament length L is on the order of the intrinsic length scale £ u . Since we have y ~ e£ u , 
we expect Eq. (2) to scale as F ~ Ae 2 /£^ ~ e 2 (A^o;) 1 / 2 , where £ is the typical value of the drag coefficient. As a 
difference, the total work done by the actuator against the viscous fluid in the period T is given by 
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which scales as W ~ e 2 T^ UJ uo 2 ~ e 2 T(A 3 ^ 5 ) 1 / 4 . To within logarithmic terms (arising in £, see below), the propulsive 
force scales therefore with frequency as F cx uj 1 ^ 2 whereas the total work done scales with a higher power as W cx cj 5 / 4 . 
For a given period and a given amount of energy available, the maximum propulsive force will therefore be obtained 
by actuating the filament only at the minimum frequency possible, uo = 27r/T, with no harmonics. 



C. Swimming 



We now consider the swimming kinematics and derive here the expected scaling for the mean swimming velocity. 
The detailed calculations will be presented in the main section of the paper. We assume that L ~ i u as this is the 
case where the optimal propulsion is expected to be generated, and we fix the amplitude, e, and frequency, cj, of the 
actuation. 

First, we consider the case where the body, of typical size a, is much larger than the filament length, i.e. a^> i^. In 
that case, the large viscous resistance of the body results in a slow-moving swimmer, and therefore small perturbations 
to the shape of filament due to swimming-induced additional drag. In that case, everywhere along the filament bending 
forces balance the drag due to the actuation and the bending amplitude, y, is on the order of e£ u . The swimming 
velocity, ?7, is then found by balancing the drag on the swimmer body by the typical filament propulsive force [15, 35] 



fiaU 



■/ 

Jo 



dy d 4 y 
dx dx 4 



and therefore 



U ~ e u 



dx ~ A 



a \og(t w /r) 



(4) 



(5) 



where the (slow) logarithmic dependence arises from the drag coefficient, £ ~ /i/ log(^,/r), where /i is the fluid 
viscosity. 

In the case of small body size, a <C t u , the expected scaling is more difficult to derive and requires a proper look 
at the local and global force balance for the swimmer. Let us denote by V the transverse velocity and by £1 the 
out-of-plane rotation rate of the swimmer, as measured as the center of mass of its body. Both V and Q average to 
zero over one actuation period but play a significant role nonetheless. Let us denote by i the typical length scale 
along the filament where bending of the filament, of amplitude y, is concentrated. The transverse force balance and 
out-of-plane torque balance on the swimmer body lead to the scalings 



lia(V + aft) - A^, iia 2 (V + aft) ~ A 



y_ 

I 2 " 



and therefore 



V - aft, 



1 



a 2 P a 3 £ 2 



(6) 



(7) 



Close to the hinge point, the drag due to the actuation is small and therefore the local force balance along the filament 
is between bending and drag due to the solid-body motion with velocity V, so that we expect 



(8) 



and therefore, given Eq. (7), 



which has the solution 



1 f I 1 



I 2 log(C/r) \ai a 2 



log 



nl/2 



(9) 



(10) 
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The deformations of the filament are therefore concentrated in a (small) region of the size of the body (to within 
logarithmic terms). The bending amplitude y, is finally determined by local force balance along the filament beyond 
the small region of size where only the drag forces balance so that Q ~ ecj, and therefore, given Eqs. (7) and (10), 
we get 

2 

(ii) 

In that case, the rotation rate exactly counterbalances the actuation and results in a straight unperturbed filament in 
the laboratory frame beyond the actuation region of size ~ £ near the actuation point. Finally, the swimming speed 
is found by balancing the drag on the filament by the propulsive force 

as the slope of the (almost straight) filament is set by the actuation (dy/dx ~ e) and the bending forces are non zero 
only in the region of size £ (with d 4 y/dx 4 ~ y 4 /£ 4 ). Consequently, we get the scaling 

1/2 

(13) 

Note that this result implies that an elastic swimmer with no head cannot swim. The results of Eqs. (5) and (13) 
show therefore that the swimmer velocity is small for both large and small body size, and therefore an optimal body 
size should exist. This will be confirmed by calculations presented in §IV. 




U ~ e uj-— 



III. SOLVING FOR THE KINEMATICS OF SWIMMING 



We now turn to the analytical calculations of the swimming kinematics of the elastic swimmer. We first present 
the assumptions used in this paper in §111 A. The intrinsic formulation of the equations of motion is derived in §111 B 
and simplified in §111 C using the small-slope approximation. Using swimming kinematics defined in §111 D, we can 
derive the free-swimming equations (§111 E), nondimensionalize the equations (§111 F) and solve for the shape of the 
oscillating filament together with the transverse velocities and rotation rates (§111 G). The values of the axial velocity 
and rotation rates are then calculated in §111 H and the final expressions for the laboratory- frame kinematics are given 
in §1111. Finally, the hydrodynamic efficiency of the swimmer is calculated in §111 J. 



A. Assumptions 

The calculations presented in this paper will be made under several simplifying assumptions. 

First, the hydrodynamics is simplified to the level of resistive- force theory [3-5, 34], a version of the equations of 
slender-body hydrodynamics [36-41] where only the leading term in an expansion of hydrodynamic forces and moments 
in powers of l/log(L/r) is conserved. In that case, the filament hydrodynamics is completely described by two drag 
coefficients, £_l and £y , relating linearly the drag forces per unit length of the filament to the local velocity relative 
to the fluid, for motion perpendicular and parallel to the filament respectively. This widely used approximation is 
asymptotically valid in the limit of very slender filaments L ^> r, and results in theoretical predictions in quantitative 
agreement with propulsive force measurements for actuated filaments [35]. To simplify the analysis, we will also ignore 
hydrodynamic interactions between the swimmer body and the oscillating filament. 

In the case of planar actuation, we then make the assumption in this paper that the amplitude of the actuation 
is small. More precisely, we will denote by e the amplitude of the oscillation of the filament slope and will derive 
the swimming kinematics in the limit where e <C 1. This has the advantage that the entire problem can be solved 
analytically and therefore shows directly the variations of swimming speeds and rotation rates with the various 
parameters of the problem. As was shown in the previous experimental study of the propulsive force mechanism by 
Yu et al. [35] , the small slope approximation gives results which agree quantitatively with numerical and experimental 
results even for large slopes, and therefore we expect the results of this paper to remain valid up to e ~ 1. 
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In the case where the actuation of the elastic filament is three-dimensional, we will furthermore neglect twist strains 
which will be generated along the filament. Obviously, the precise nature of the twist strain density in the filament 
depends on the exact actuation mechanism at the base of the filament. However, it is possible to show that such twist 
strains do not influence the kinematics of the swimmers in the limit which is of interest to us. Indeed, as a difference 
with bending amplitude which hyperdiffuses as described in Eq. (1), twist strains satisfy a diffusion equation, and 
the intrinsic twisting- hydrodynamics length scale is = (C/o;^) 1 / 2 , where C is the twisting modulus and £ r is the 
rotational drag coefficient for the filament, that is £ r = Air/nr 2 [28]. Since the relevant filament length for effective 
swimming is L ~ and since for most solids A ~ C, we have L/L^ ~ l^/L^ ~ (r/£ a; )(log(£ a ;/r)) 1 / 2 . As the 
filament is slender, the typical twist penetration length scale is therefore much larger than the relevant filament 
length, L/L^ <C 1, and as a consequence twist strains are always in diffusive equilibrium: They vary linearly between 
zero at the free-end of the filament and a constant at the base, which can be obtained by a local balance between 
viscous torque and twist. Now, with this balance, it is possible to show [28] that the twist term in the equation of 
motion for the filament is smaller by a factor r 2 /£ 2 (within logarithmic terms) than the bending term, and therefore 
can safely be neglected. This also means that we can ignore in this paper possible buckling (whirling) instabilities, 
which occur above a critical rotation frequency u c ~ A/( r L 2 ~ A/fir 2 L 2 [28]. Indeed, in the case where L ~ we 
get ijj c jijj ~ £ 2 /(r 2 log^/r)) ^> 1, and therefore buckling instabilities appear for much larger values of the typical 
actuation frequency and can be neglected. In that case, we will also ignore local rotational drag along the long filament 
[27]. 



B. Intrinsic formulation and equations 



The elastic energy of the flexible filament as a function of its confirmation is given by 



2 Jo 



A 



-JO" 



ds 



(14) 



where A is the bending modulus, a(s) the Lagrange multiplier (tension) which enforces local inextensibility and s 
the curvilinear coordinate along the filament (0 < s < L). As discussed above, twist strains are not included in this 
paper and therefore do not appear in Eq. (14). In the case of filament of circular cross-section with radius r, we 
have A = 7rr 4 £^/4, where E is the Young's modulus of the material composing the filament. Assuming resistive force 
theory for the fluids forces (see above), the intrinsic formulation for the filament elasto- hydrodynamics is obtained by 
calculus of variation from Eq. (14) and is given by 



G« + -«>!- = "4? + s('S). 



(15) 



where u is the local instantaneous velocity along the filament. The equation for the Lagrange multiplier <r, which 
enforces inextensibility (r s • r s ) t = is given in an implicit form by 



du 

These two equations have boundary conditions which are given by 



OS 6 OS 



Q2 r 

xt + i— 



Q3 



dr 



Fext = -A-;-; +CT — , tX 



ds 3 



x t - A 



ds 2 
d 2 r 



= 0, at s = 0, 



ds 2 



0, at s = L, 



(16) 

(17a) 
(17b) 



where F ext and T ext are externally applied forces and torques at the ends of the filament and t = r s 



C. Small-slope approximation 

Let us consider a cartesian coordinate system {e x , e y , e z } moving with the swimmer and located at the base of the 
elastic filament such that e x is directed in the mean direction of the filament (see Fig. 1). The actuation oscillates 
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therefore around the e x axis. The laboratory coordinate system is defined as {ei,e2,es} and is chosen such that the 
average swimming occurs in the ei direction. The filament positions is described by the functions y(x,t) and z(x,t). 
We assume that the filament shape is slowly varying, that is 



dy 
dx 



«1, 



dz 
dx 



«1, 



(18) 



so that a point on the filament is defined as r « xe x + r_L, where = y(x, t)e y + z(x, i)e z , and we also have s « x. 
In that case, Eqs. (15) and (16) become, when written in the frame moving with the swimmer, 



. f . . .dy , , . v d z da 



dy 

(£|| ~ €-Oa- u x + Z,±Uy + (f y - £_l) 



9x 

dz dy dz 



dx dx 



, <9 4 v 9 
[ — - -| 

dx 4 dx 
d 4 z d 



dx 2 



- 1 



d 2 r± d 4 r± A^ T± ^ 5r ^ 



dx 2 dx 4 



dx dx 5 





(19a) 


dx) 


(19b) 




(19c) 
(19d) 


dx) 



where, {u x , u yi u z } are the local velocity components of the filament. The leading-order boundary conditions at either 
end of the filament are given, for all times, by 



^d 3 r± dr± 

dx 3 dx 
. d 3 r± 



cre x 
dr ± 



Te ^ y ~ A dx 2 ' 



Text, z 
d^ r 

L dx 2 ' 



{ d 2 y 

dx 2 ' 



at 5 = 0, 



Fext = - A ^ ey + a ^x- +ae ^ T ^ y = ~ A d^ 
The most general periodic actuation at the frequency co is 

2/(0, t) = 0, ||(0,t) =ecoaojt, 
dz 

z(0,t) = 0, — (0, t) = <5esint^. 



t ext, z 



= A 



d 2 y 
dx 2 ' 



at s = L. 



(20a) 
(20b) 

(21a) 
(21b) 



where < 6 < 1 measures the extent of the three-dimensionality of the actuation. In that case, 6 = represents a 
purely planar actuation, whereas for 6 = 1, the actuating hinge sweeps a cone. With these notations, the small slope 
approximation of Eq. (18) is written e«l. 



D. Swimming kinematics 

Let us denote by U(t) the instantaneous velocity of the swimmer body and fl(t) its instantaneous rotation rate 
around the hinge point [43] . The motion of the filament is then described as a superposition of a solid body translation 
at velocity U, a solid body rotation with rotation rate f2, and a relative motion due to the oscillations of the filament 
and described by the functions dy/dt and dz/dt. As a consequence, the local velocity components along the filament 
are given by 

u x = U x + zQ y - yQ z , (22a) 

u y = Uy + ^ + xft z - zQ x , (22b) 
dz 

u z = U z + — +yfl x - xtty- (22c) 



E. Free-swimming assumption 



We consider a free-swimmer in Stokes flow, so the total force and torque on the swimmer must vanish. This means 
that the forces and torques in Eqs. (20a) have to balance the hydrodynamic forces and torques on the body of the 
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swimmer, whereas the forces and torques in Eqs. (20b) must vanish. In this paper the swimmer will be assumed to 
have an axisymmetric body around the x axis. This allows the mathematical formulation to remain manageable and 
the generalization to more complex body shapes is straightforward. In the center of mass of the swimmer body, the 
resistance matrices for fluid forces and torques acting on the body are written, by symmetry, as 



Fpyf — — 





(23) 



where everything is measured relative to the center of mass of the body. In the frame of reference located at the base 
of the filament (see Fig. 1), these relations are now written as 



(24) 




ft, 



(25) 



a 2 TZi u 



where a is the distance between the center of mass of the swimmer body and the actuation point (see Fig. 1) and 
ft the rotation rate around the actuation point. Using these notations, the force-free and torque-free conditions at 
s = 0, Eq. (20a), are written as 



T$ u U x {t) = "(0,*), 



nl u u y {t)-anl u n z {t) 

-aTli u U y {t) + (TZi Q + a 2 Tli u )n z {t) 

ani u U z {t) + (TZ^ 1 + a 2 1li u )Vl y {t) 
while the conditions at s = L (Eq. 20b) become 



4?(o, t) , 

d^ z dz 
d 2 z 



dx 2 



(L,t) 



dx 2 



(L,t) 



dx 3 



(L,t) 



d 3 z 
dx 3 



(L,t) = a(L,t) = 0. 



(26a) 
(26b) 
(26c) 
(26d) 
(26e) 

(27) 



The system of partial differential equations is closed by writing down the overall torque balance in the ^-direction, 
leading to an equation for Q x 



f 

Jo 



Aly 



d A z 



d 4 y 



dx 4 dx 4 



d 



dy 



dx V dx 



y- 



d 



dz\ 



dx V dx J 



dx = A 



dy d 2 z dz d 2 y 
dx dx 2 dx dx 2 



(28) 



x=0 



where we have used integration by parts and conditions (21) and (27). We have now as many equations as we have 
unknowns. There are nine unknowns, {U x ,Uy,U z ,Q x ,Q y ,£} z , cr,y, z} and nine equations (Eqs. 19a, 19b, 19c, 19d, 
26b, 26c, 26d, 26e, 28). The equation for a is second order, and is accompanied by two boundary conditions (Eqs. 26a 
and 27). The equations for y and z are fourth-order and are accompanied by four boundary conditions each (Eqs. 21 
and 27). 



F. Nondimensionalization and simplifications 

We now nondimensionalize the equations of motion. We scale lengths by the intrinsic length scale time by 
a; -1 , rotation rates by a;, velocities by a;^, resistivities by £±£™ (n = 1 or 3 depending if it is a force- velocity or a 
torque-rotation rate resistivity), forces by ^±£^uj and torques by £±£^uj. 
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Furthermore, we can use the fact that e is small to simplify the dimensionless equations further. Let us evaluate 
the leading order power of e for each of our nine unknowns. From Eq. (21) we see that y ~ z ~ e. Consequently, from 
Eqs. (26b), (26c), (26d), and (26e) we see that U y ~ U z ~ Q y ~ Q z ~ e. We then get from Eq. (28) that Q x - e 2 and 
Eq. (19d) shows that a ~ e 2 , so that Eq. (26a) leads to U x ~ e 2 . The magnitude of the axial swimming speed and 
rotation rate are therefore one order of magnitude smaller than the transverse velocities and rotation rates. These 
scalings allow one to simplify the equations for y and z further, and we obtain, using the same symbols as for the 
dimensional variables for simplification, 



-i + U y (t)+xQ z (t) = (29a) 



dv d^v 

-JL + Uy{t ) + X n z (t) = 

37 3^ 7 

g- t +U z (t)-xn y (t) = (29b) 



Because of the e scaling for the transverse problem is an order of magnitude larger than the scaling e 2 for the axial 
problem, we see that the axial unknowns (U x ,£l x ) have disappeared from Eq. (29). We can therefore solve these two 
problems in two separate times. First, we solve Eq. (29) for the filament shape {y, z} and the transverse swimming 
kinematics {U y ,U z ,Q y ,£l z }, with boundary conditions 

2/(0,*) = 0, ||(0,t) = ecost, = 0, = 0, (30a) 

3z 3 2 z 3 3 z 

*(0,t) = 0, -{Q,t)=8csmt, —(L,t)=0, —(L,t) = 0, (30b) 

where L now refers to the dimensionless length of the filament. The dimensionless resistance equations become, with 
all the symbols referring now to dimensionless variables, 

ni u U v {t)-aRi u n z {t) = "0(0,0, (31a) 

-ani u U y (t) + (ni Q +a 2 n^ u )n z (t) = f|(0,0, (31b) 



nl u U z {t)+aRl u ^l v {t) = "0(0,0, (31c) 
a1li u U z {t) + {Tlf + a 2 1li u )n y {t) = "0(0,0, (31d) 



which can be inverted to give 



^) = -(^^)S(0.0 + ^S(0, f ). (32a) 



= - ( ni L + ut?F ) S(0,0 + iS(0,0, (32c) 



(32d) 

Once we have the solution for {y, z, U y , U z , Q y , Q z }, we can use Eqs. (19a), (26a) and (28) to get U x and Q x . Note 
that the dimensionless version of Eq. (19a) is written as 

/. x / dy dz \ da ^ ^ 

«* + (W)(^ + ^)=7^ (33) 

where 7^ = Integration of Eq. (33) along the filament, using Eqs. (26a) and (29), leads to the formula we will 

use for the axial swimming velocity, U x , as 



(7- 



_K( u + L)U x = n z / ydx-Q y / zdx^(l-j ± ) 
Jo Jo 



1 f d 2 y\ dy d 3 y 1 fd 2 z\ dz d 3 z 

2 \ dx 2 J dx dx 3 2 \ dx 2 J dx dx 3 



■ (34) 

X = 



10 



It is important to note at this point that the small-slope approximation has resulted in a partial simplification of 
the problem: The axial velocity and rotation rate of the swimmer being one order of magnitude smaller than the 
transverse velocities and rotation rates, the axial swimming kinematics is slaved to the transverse kinematics. The 
problem of determining the filament shape and the transverse swimming kinematics cannot, however, be simplified 
any further and both still have to be solved simultaneously. 



G. Solving the transverse problem: Filament shape and swimming kinematics 



Let us now solve Eq. (29). Since the forcing is harmonic, we will solve these equations in Fourier space and write, 
for all variables, A(x,i) = ?R{A(x) exp(— it)}. Using Eq. (32), we have the relations 

m)+ «. m , _ ( »yy ) g ( o,o + <-£>0<o, ( >, 



U Z (t) -XSly(t) 



+ a(a + x)K I [ u \ d^z (a + x) d 2 z 

n™iV; a J dx* K ' ' + dx* { ' '■ 



(35b) 



From Eqs. (29) and (35) we see that the generic equation satisfied by both y and z is a hyperdiffusion equation 
forced by a first order polynomial whose coefficients depend on the boundary condition of the solution. This integro- 
differential equation is a consequence of the problem-coupling discussed in the introduction of the paper and reflects 
the nonlocal aspect of locomotion without inertia where, at all times, velocities and rotation rates adjust so that total 
forces and torques sum up to zero. 

Let us denote by Q{x; (3, a, A, ji, h), the solution to the differential equation 



d 4 1 d 3 C d 2 C 

- i+ Ax~A C(X) = + ax) d^ (0) + (A + ^ } d^ (0) ' 



(36) 



with boundary conditions 



C(o) = o, fjo) 



d^c 

dar 3 



(h) = 0. 



In that case, it is easy to see from Eqs. (29) and (30) that, if we define 

_ / 1 a 2 a a 1 



then we have 



(37) 



(38) 



y(x,t) = e^ie-'X^x)}, z{x,t) = Se^ie^Uix)}, 



(39) 



and we note that z = iSy. 

The analytical solution to Eqs. (36)- (37) is given by 



C(x) = ^ A n e a - X + Bx + C, 

n=0 



where a n = exp(i(l + 4n)7r/8) (0 < n < 3) and where the six constants satisfy the linear system 



(40) 



1 

a 



ale aoh 



ai 



a-i e 



1 10 1 

1 

ale aih a\e^ h a\e^ h 

aKX + fiao) a\(\ + f3ai) ^(A + ^o^) ol\{\ + Pas) i 
aKfi + aao) a\{n + aai) aKfi + cra 2 ) + eras) i 



OL2 

a\e^ h 
a\e^ h 



a 3 
a\e^ h 
a\e^ h 





-A ~ 




"0" 




A x 




1 




A 2 









As 









B 









C 








(41) 
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With this solution known, we get the transverse velocities and rotation rates 



^ = ~ e C "' (0) + # C ' (0) ' (42a) 

= -^rCr(0) + ^r^(0), (42b) 

U z = iSU yi (42c) 
Cl y = -iSCl z . (42d) 

This completes the solution of the transverse problem, and we now have the expressions for both the filament shape 
(Eq. 39) and the transverse swimming kinematics (Eq. 42). 

H. Solving the axial problem: Swimming velocity and rotation rate 

We can now solve the second problem, namely find the expression for the axial velocity and rotation rates. Using 
the Fourier-space notation defined above, it is then easy to see that the axial swimming velocity (Eq. 34) has two 
components, U x = (U x ) + U' x , where (U x ) is the steady component and U' x the zero- mean oscillatory component, 
which are given by 

(Ux) = 2 ^tfl L) » { m l + (1 - 7±) (^ICI(O)I 2 ~ Ci(0)C (0)*) } , (43a) 



m z + (i - 7± ) Qcl(o) 2 - CKo)cr(o))] } , (43b) 



where we have defined X = J Q (±(x) dx and Q z = Cl z /e (which is of order one). Note that when 5 = 1, the unsteady 
component is exactly equal to zero (although the trajectory is not - it is a helix). This is because in this case, the 
shape of the filament is steady in the frame rotating with the hinge [27]. 

The equations above can be simplified further by noting that the integration in space of the Fourier-transform of 
Eq. (29a), using Eq. (31a) as boundary condition, leads to the value of X as given by 

yFU ^ \ -(-nFU 



X = iQ z \aKi u - — j - i(Ki u + L)U y , (44) 
where U y = U y /e. As a consequence, we get 

®{m* z } = {n F A u + L)z{u y ti* z } (45) 

and, given Eq. (42), it is straightforward to show that 

= n Fu n Ln Z{Cl(0)Cf(0y}, (46) 

so that the expression for the axial swimming velocity is simplified to 
e 2 (l + <5 2 ) 



(U x ) 



(47a) 



(i - 7± )5R |i|c!(o)| 2 - Ci(0)cr(0)*} + ( ^I^lo ) s{C"(o)C(o)*} 

K = 2{ ^fV 2 +L) [ (1 " 7±)5R { e ~ m (^ C " (0)2 " } < 47b ) 

+(nl u + l)9 {u v n z e- 2U } + ( y " an ± u ) 9 {^ e " 2it } • 

Finally, the leading order version of torque balance in the x direction, Eq. (28), shows that Vt x does not have any 
oscillating components (i.e. is steady) and is equal to 

^ = (fix) = ^3{Ci(0)Cr(0)}, (48) 

which finishes the solution for the axial problem. We now have the expression of all three-components of swimming 
velocities and rotation rates in the frame moving with the swimmer body. 
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I. Swimming kinematics in the laboratory frame 



The final step in the calculation is to transform the expressions we have for the swimming kinematics in the body- 
fixed frame to the laboratory-fixed frame. The general calculation for this transformation is presented in Appendix A. 
It is worth noting that the distinction between the two frames of references has rarely been discussed in the context 
of swimming micro-organisms but can have important consequences if not properly taken into account. In our case, 
using the transformation given by Eqs. (A8) and (A9) as well as the swimming velocities obtained above, we find the 
final formulae for the leading-order swimming speeds in the laboratory frame as given by 

■R' r , 

/-ft I 



(Ui) 

(U2) 

(Us) 



U' 2 

u 3 



e 2 {l+5 2 ) 
2( 7 ±ftf + £) 
0, 
0, 

e 2 (l - 5 2 ) 
2(^11™ + L) [ 



(i-7j5R{^IC"(o)| 2 -Ci(o)C"'(or 



+ 



(o)cr(o)*} 



(l-7±)» 



<!(o) 2 -ci(o)cr(o) 



)} 



+{Ki U + 1±TZ^ U + 2L)3 {U y n z e- 2it } + (y - a7 ^) 9 & 



-lit 1 



[cos(n x t)$l {e- u U y } 
[sin({l x t)*t{e- it Uy} ■ 



- Ssm(n x t)^l {e'^Uy}] , 
5 cos(n x t)3 {e- u U y }] . 



,(49a) 

(49b) 
(49c) 

(49d) 



(49e) 
(49f) 



The expressions given in Eq. (49) are the main results of this paper. Note that when the drag is isotropic (7^ = 1, 
1Z FU = TZ^_ U ), the mean swimming velocity is exactly equal to zero. Drag anisotropy is therefore crucial for locomotion 
without inertia [9]. Note also that in Eqs. (49e) and (49f), time appears as a parameter as a result of the separation 
of time scales in the limit of small e (see discussion in Appendix A). 



J. Hydrodynamic efficiency 



We define in this paper the efficiency of the motion, £ , as the ratio of useful work (defined as the work necessary 
to move the entire swimmer at the steady speed (Ui)) by the total work done by the swimmer, 



<Wi> 



da; 



which becomes in dimensionless variables and at leading order in e 2 



£ = 



nr + U) {Ul) 



/ 

JO 



dx 



The first term in the denominator is given by 

u f + n l = nl u (u y - an z ) 2 + n F A u {u z + an y ) 2 + izi n n 2 y + ni n n 2 z . 

while the second term can be evaluated in Fourier space and we obtain the efficiency as given by 



2 n 



FU 



€ 2 (l + d 2 ) 



n I [ u \u y -an z \ 2 + n I [ n \n z \ 2 + I |c""(^)l 2 d^ 

Jo 



(50) 



(51) 



(52) 



(53) 



with (Ui) given by Eq. (49a). Note that since (Ui) scales linearly with (1 + S 2 ) and appears squared in Eq. (53), we 
obtain the result that a three-dimensional actuation (S 7^ 0) is always more efficient than a planar one. 
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Velocity (I) Velocity (II) Velocity (III) Efficiency (IV) 
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L 


\(Ui)\ 




Wi)\ 


S 




L 




e 2 


e 2 (2a + L) 


2e 2 a 


e 2 








(xlO 3 ) 


(xlO 3 ) 


(xlO 3 ) 




Spherical body 














Optimal I 


n on 


Z.Ol 


o.y 


L.I 


4.0 


U.izz/o 


Optimal II 


0.29 


2.38 


6.7 


1.8 


4.9 


r\ -\ r\r\Cr/ 

0.100% 


Optimal III 


0.18 


2.83 


5.7 


1.5 


5.7 


r\ r\ r-7f \ ()-/ 

0.079% 


Optimal IV 


0.37 


2.70 


6.1 


1.3 


3.1 


0.142% 


Elongated body 














Optimal I 


0.62 


2.78 


17.9 


2.9 


5.2 


0.418% 


Optimal II 


0.44 


2.82 


16.7 


3.1 


6.7 


0.369% 


Optimal III 


0.24 


3.34 


12.6 


2.6 


7.9 


0.336% 


Optimal IV 


0.49 


3.20 


16.9 


2.6 


5.4 


0.465% 



TABLE I: Geometrical and actuation characteristics of the optimal elastic swimmers (spherical body and elongated body of 
aspect ratio 500). Four different quantities are optimized: The velocity in the laboratory frame (I), the velocity in swimmer 
length per unit beat (II), the velocity in body length per unit beat (III), and the mechanical efficiency of the swimmer (IV). 
Velocities and efficiencies are given for the planar swimmer (8 = 0) but the geometrical characteristics of the optimal swimmers 
are independent of the value of 8. The swimmers are illustrated in Fig. 2. 



K. Asymptotic limit of long filament 

With our analytical formulae, we can now derive the swimming velocity in the limit of a long filament L ^> 1 
(that is, L ^> £ u in dimensional variables). Although we expect the swimming velocity to decrease to zero in this 
case, it is the biologically relevant limit for the motion of spermatozoa [5, 15]. Obviously, spermatozoa use a different 
swimming mechanism as the one described in this paper, so our purpose is merely to be able to compare swimming 
performances. 

In the limit of large body a ^> 1 (a ^> in dimensional variables), it is easy to see that the mean velocity, given 
by Eq. (49a), becomes 

™-^(^). 

In the limit of small body a <C 1 (a <C ^ m dimensional variables), one needs to write down Taylor expansions for 
each of the body resistivities in cl/£ (jJ ^ which is tedious but straightforward, and the mean swimming velocity is given 

by 



IV. PERFORMANCE AND OPTIMAL DESIGN OF THE ELASTIC SWIMMER 



As we have now analytical solutions for the complete swimming kinematics, we can study the performance of the 
elastic swimmer. 
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III 
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FIG. 2: Optimal swimmers. Left (blue): Optimal swimmers with spherical body. Right (red): Optimal swimmers with 
elongated body, (see caption of Table I for description). 



A. Optimal spherical swimmers 

As we have noted before, and as is confirmed by our analytical formulae, swimming velocities and efficiencies 
decrease to zero for both large and small filament length and body size. As a consequence, optimal swimmers exist. 
The geometrical characteristics of these optimal swimmers are given in Table I and represented schematically in Fig. 2. 
For the calculations presented below, we have assumed the filament to be far from any boundaries and therefore the 
appropriate drag coefficients are given by 

log(L/r)-l/2' U log(L/r) + l/2' 1 } 

and we have chosen the aspect ratio of the filament to be L/2r = 500 [44]. We will furthermore assume the body to 
be a spheroid of revolution [42]. 

We first determine the optimal swimmers with spherical body (Table I and Fig. 2 left). There are a variety of 
ways to define the cost function to optimize, and we study four different optimality conditions. The first one is the 
swimming speed (Ui) (I). The second measure of performance is the swimming speed expressed in swimmer length per 
beat, that is (U\)/(2a + L) (II) . The third important velocity is the swimming velocity expressed in body length per 
unit beat, and is given by (Ui)/2a (III). Finally, we also consider the swimming efficiency as a measure of performance 
to optimize (IV). As shown in Table I and Fig. 2, the geometrical characteristics of the optimal swimmers are a strong 
function of the performance index which is chosen (but are independent of the value of S). In all cases, the optimal 
filament length L is on the order of the intrinsic length scale i^. 

The variations of the swimming speed and efficiency of the optimal swimmer I (the fastest) and IV (the most 
efficient) with changing filament length and body size are displayed in Fig. 3. In particular, we see that the swimmer 
cannot move if either its filament or its body is too small. In the limit of large swimmer, the velocity also decreases 
to zero as the inverse of the swimmer size. 

It is interesting to note that, out of the three length scales which can be a priori chosen independently in designing 
an elastic swimmer - namely a, L, and - their relative magnitude is fixed for the optimal swimmers (see Table I) 
and therefore only one of them can be chosen arbitrarily. For a given performance index, once a specific body size or 
filament length or actuation frequency (through I J) is chosen, everything else is fixed and there is only one possible 
optimal swimmer. 
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FIG. 3: Variation of the mean swimming velocity, (Ui), and the mechanical efficiency, £, as a function of the body to filament 
length ratio, a/L, and the dimensionless filament length, L/£ u , for the optimal swimmer with spherical body I and IV (see 
Table I). The values are displayed for the planar swimmer (6 — 0). Top: Variation of the swimming velocity for the fastest 
swimmer (a/L — 0.3, L/t^ — 2.51); Bottom: Variation of the efficiency for the most efficient swimmer (a/L — 0.37, L/t^ — 2.7). 



B. Optimal swimmers 



We then study how the shape of the swimmer body influences the swimming performance and we find that a better 
performance is always obtained for a long slender body of large aspect ratio in the swimming direction (the overall 
best is obtained in the limit of an infinite aspect ratio). We present in Table I the performance of the optimal swimmer 
with an elongated body of aspect ratio 500, same aspect ratio as the filament. Both swimming speed and efficiency 
improve significantly by taking a slender swimming body. The optimal swimmers are displayed in Fig. 2 (right). 
Again, the choice of a performance index has consequences on the resulting optimal shape. 



C. Driving torque 



From a possible practical standpoint, it is important to quantify the internal torque necessary to drive the actuation 
at the given amplitude, e, and frequency, uo. We consider here the case of planar actuation. In that case, the oscillating 
torque is equal to the torque given in Eq. (26c) and therefore, in dimensionless form, the torque amplitude is given 
by \T\ = e|C"(0)|. The variation of the internal torque with the swimmer size is illustrated in Fig. 4 in the case of the 
fastest swimmer with a spherical body (optimal swimmer I). 



D. Filament shape 

Finally, we can use our analytical solution to characterize the shape of the flexible filament as the swimmers moves, 
and compare it to the shape obtained when the filament is actuated but does not swim [15, 35]. The results are shown 
in Fig. 5 for the planar swimmer with a/L = 0.3. Since the swimmer body can always rotate to relieve some of the 
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a/L L/iu 

FIG. 4: Variation of the amplitude of the internal torque supplied by the actuating mechanism as a function of the body to 
filament length ratio, a/L, and the dimensionless filament length, L/£ u , for the fastest swimmer with spherical body (a/L — 0.3, 
L/iuj — 2.51, see Table I) and for planar actuation (5 = 0). 



applied torque at the base of the filament, the curvature of the filament is smaller in the free-swimming case than 
it is in the case where the filament is not free to move (compare the shapes on the left and at the center of Fig. 5). 
We also note that the filaments in both cases display exponentially decaying amplitude, a feature which makes this 
problem intrinsically different from eukaryotic flagellar propulsion where large-amplitude oscillations are present along 
the entire filament. 



V. ELASTIC SWIMMING WITH MORE THAN ONE FILAMENT 

In this final section, we discuss possible improvements on the design of elastic swimmers. One drawback of using 
a single filament as a propeller is the oscillatory nature of the swimming kinematics. As a result, a lot of effort goes 
into propelling the swimmer body in a direction which is different from the main swimming direction. In fact, in the 
limit of small actuation amplitude e, we have seen above that the mean swimming speed is of order e 2 whereas the 
transverse swimming speed (and rotation rate) are both of order e and therefore one order of magnitude larger. We 
propose in this section to use more than one filament in order to have better control over the instantaneous swimming 
direction. For simplicity, we will ignore hydrodynamic interactions between the filaments in the analysis below. 



A. Two filaments in planar motion 

If the swimmer has two filaments which are (a) positioned exactly symmetrically with respect to the axis of symmetry 
of the swimmer body and (b) are actuated with 180° out of phase, then all transverse forcing due to one filament 
will be exactly canceled by the second one, and this will result in a straight (yet unsteady) swimming trajectory (see 
Fig. 6a). From a mathematical standpoint, the dynamic in the x and y directions will be coupled for each filament, 
and we expect therefore a steady swimming of order ~ e 2 with oscillations of magnitude ~ e along the same direction. 
In that case, the only nonzero component of the velocity is U(t), along the ei direction (the direction e x still refers 
to the direction of each of the filaments) and there is no rotation rate. The equation for the shape of each filament 
becomes 



dy 
dt 



dx A 



UsinO, 



(57) 



where is the angle between the average position of the filament base and the swimming direction. Force balance in 
the x direction gives 



(7±W| + 2L cos 2 0)U(t) =27_Lsin<9 



d 3 y 



dx 3 



2(1 -7_l)cos0 



x=0 



1 (<Py 

2 \dx 2 



dy d 3 y 
dx dx 3 



(58) 
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FIG. 5: Shapes of the elastic filament as a function of the dimensionless length L/iu. Left (black): Shape of the filament 
in the case where the actuation point is fixed in space. In that case, the filament does not swim and its shape is that found 
theoretically by Wiggins and Goldstein [15] and experimentally by Yu et al. [35] . Center (blue): Shape of the filament in the 
case of the free-swimming as seen in the frame {e x ,e yi e z } translating and rotating with the swimmer body. The body size is 
a/L = 0.3. Right (red): Shape of the filament for the same swimmer as seen in the laboratory frame {ei,e2,es}. In that case, 
the instantaneous swimming velocities have been subtracted in order to be able to compare shapes. 



This is solved order by order for both y(x, t) = eyi + e 2 i/2 + ... and U{t) = eJJ\ + e 2 /72 
time average. At leading order we find that y\ is solution to the equation 



dyi 
dt 



27 ± sin 2 



d 4 yi = 
dx A -f ± 7Z^ u ^2L cos 2 ( 



d 3 yi 



dx 3 



with boundary conditions 



2/1 (0, *) =0, 



9yi x d 2 yi 
— (0, t) = cost, 



dx 



dx 2 



(L,t)=0, 



J x=0 



d 3 yi 
dx 3 



and only U2 has a non-zero 
(59) 



(L,t) = 0. 



(60) 
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(a) (b) 

FIG. 6: Swimming with more than one elastic filament; (a) Swimming with two symmetric filaments in planar motion leads 
to unsteady straight swimming (§V A); (b) Swimming with three symmetric filament pairs (so six filaments) in planar motion 
with each 2tt/3 phase difference leads to steady straight swimming (§VB). 



Consequently, y\ is given by 

y 1 (x,t)=$t{e- it ( 1 (x)}, Ci(*) = CU 



27 ± sin 2 



7± ftf 7 + 2Lcos 2 6> 

and then the leading order component of the velocity can be calculated by 

2~/± sin 



,0,0,0,L 



Ui(t) 



J±n£ u -f 2Lcos 2 <9 



5R{e-^C(0)}, 



(61) 



(62) 



which is purely oscillating (no mean component). We need to go to next order to calculate the mean component of 
the swimming. Since there is a component of the mean swimming velocity normal to the filament, each filament is 
therefore slightly asymmetric (the asymmetry, i.e. the mean value for yi, is of order e 2 ), a feature which is not present 
when we consider the case of a single filament. The second order term, y^ satisfies the equation 



dm d 4 y 2 
dt dx A 



U 2 (t) sin (9, 



(63) 



with homogeneous boundary conditions. Since U2 is expected to have both a mean value and an oscillating part, 
U2(t) = (U2) + U' 2l we consider them separately. The oscillating part of U2 will lead to an oscillating value for y2 
which will average out to zero in Eq. (58), and therefore we do not need to solve for it. On the contrary the steady 
component of U2 leads to a steady deflection of the filament, described by 



d 4 (t/ 2 ) 
dx 4 



= ([/ 2 }sin6>, 



with homogeneous boundary conditions. This can be integrated to give 

„4 ,,,3 



{y 2 ){x) = {U 2 )smOf{x), f(x) 



24 



x~L 
6~ 



x 2 L 2 



(64) 



(65) 



which can be then used to evaluate the mean swimming speed by averaging Eq. (58) at second order, and we find 



(7-L^-il +2Lcos 2 6» + 2j ± Lsm z 9) (U 2 ) = (1 -7_l)cos6» 



-|Cf(0)| 2 -5K{CUO)CP(0)} 



(66) 



When the swimmer does not have a body, T^^ u = 0, we obtain the "elastic scallop", flexible version of the two-arm 
swimmer discussed by Purcell and whose name is at the origin of the so-called scallop theorem [8]. Note that a 
two-filament swimmer with no head can swim, whereas a swimmer with a single filament cannot swim without a 
head. In the case where L ^> 1, the function (j can be simplified further [15] and we find the average swimming 
velocity as given by 



(U) 



y/2-1 



e 2 (l — 7_l) cost 



V2 l±nf\ u + 2L cos 2 + 2 7± L sin 2 ' 



(67) 



with corrections exponentially small in L/£ u 
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B. Steady swimming: more than two filaments 



The setup with two filaments allowed us to get motion on a straight line. However, we obtained a steady swimming 
speed of order ~ e 2 with longitudinal oscillations of order ~ e, which is not ideal. Here we ask the question, can we 
design a swimmer which moves on a straight line steadily? Since we look for straight motion, we will always look for 
pairs of filaments, located symmetrically with respect to the body symmetry axis (recall that the body is assumed to 
be axisymmetric), and mirror image to one another (see Fig. 6b). We need to find the minimum number of filament 
pairs together with the phase differences between each of them which is necessary to choose in order to achieve an 
overall steady propulsive force despite the unsteady propulsive forces (and shapes) of the individual filaments. When 
appropriately chosen, these phase differences will cancel out all of the ~ e contributions from each individual filaments 
and should leave us with a steady ~ e 2 swimming speed. 

The equation for each of the individual filaments in one pair (1 < n < TV, where N is the number of pairs) is given 

by 



with the boundary conditions, 

y n (o,t) = o, 



dy n 
dx 



dy n d A y n 
dt dx± 



(0,t) = d?{c n e-* £ }, 



E/sin0, 



d 2 y n 
dx 2 



d 3 y n 
dx 3 



(L,f) = 0, 



(68) 



(69) 



where the set of complex constants c n is determining the phase difference in the actuation for each filament pair. Note 
that the filaments could be located anywhere around the swimmer body, but the assumption of ignoring hydrodynamic 
interactions will be the most valid when they are the furthest apart, so we assume the filaments to be distributed with 
an angle tt/N apart from each other. Assuming a steady swimming speed, the solution to Eq. (68) is given formally 
by 



y n (x, t) = U sin 0f(x) + e3?{c„e u g(x)}, 



(70) 



where / is the steady deflection above (Eq. 65) and g(x) is the Wiggins & Goldstein shape (shape of the elastica when 
it is actuated at one end and does not swim), that is 



g(x) =C(z;0,0,0,0,L). 
The force balance in the swimming direction leads to the swimming speed as given by 



(71) 



(7±ftf + 2iV cos 2 9L)U = 2 7 _l sin 9 



N 

E 

n=l 



d 3 y n 
dx 3 



+ 2(1 -7 ± )cos6> 



x=0 



N 

E 

n=l 



d 2 y n 
dx 2 



dy n d 3 y n 
dx dx 3 



(72) 



J x=0 



Given that the filament shape appears in Eq. (72) with linear and quadratic terms, it is straightforward to see that 
the smallest number of filament pairs necessary to achieve straight steady swimming is such that 



N 
n=l 



o, 



N 

E c « 

n=l 



0. 



(73) 



The minimum number of filament pairs is therefore N = 3 (so 6 individual filaments) with c n — e 2m7r / 3 (n = 1, 2, 3). 
The value of these constants show that if one pair of filaments displays base oscillations of the form cost; then the 
other two pairs need to oscillate as cos(r; ±2tt/3). In that case, we obtain that the (steady) swimming speed is given 
by the quadratic equation a quadratic equation 



(-y±Kf u + 6Lcos 2 <9 + 67 ± Lsin 2 0)U = -(1 - J±)U 2 L 4 sin 2 (9 cos (9 + 3e 2 (l - 7± ) cos (9 
and since we have assumed e <C 1 we get the leading order solution given by 



i| 5 "(0)| 2 -K{ 5 '(0V"*(0)} 



(7-L^jl + 6L cos 2 9 + Qj ± L sin 2 9)U = 3e 2 (l - 7_l) cos 9 



2l/(0)| 2 -^W"*(0)} 



(74) 
(75) 
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Eq. (75) is the only instance where, up to geometric factors, the equation for the swimming speed is simply given by 
a balance between the viscous drag on the filaments and the propulsive forces from Wiggins & Goldstein for all six 
filaments. In the limit L ^> 1, the approximate solution for the swimming velocity is now given by 

3(>/2-l) e 2 (l- 7± )cosfl 

y/2 l±Kf\ U + 6Lcos 2 (9 + 67 ± Lsin 2 (9 1 } 



VI. DISCUSSION 

We have presented in this paper an analytical treatment of the locomotion of an elastic swimmer in the limit of 
small amplitude actuation. This is arguably the simplest device which takes advantage of the coupling between drag 
and bending forces for locomotion purposes. Our study is different from and improves upon previous work in many 
aspects: (a) Our analysis does not violate force balance nor torque balance and the constraints of free-swimming 
motion are fully enforced; (b) We include in our study the distinction between body-fixed and laboratory-fixed frames 
of references and calculate the swimming kinematics in both frames; (c) The coupling between the two problems - 
finding the shape of the filament and the swimming kinematics - is discussed for the first time and is solved in a 
self-consistent manner; (d) Our study produces analytical formulae; (e) We obtain the expected limits of vanishing 
swimming speed for small and large swimmer sizes; (f) We characterize the geometry and performance of optimal 
swimmers; (g) The study of elastic swimming with more than one filament is presented, as it allows a better control 
on the swimming trajectories. 

It is enlightening to compare the performance of the elastic swimmer with that of swimming micro-organisms - 
which of course use different actuation mechanisms. First, let us consider the bacterium E. coli [6]. The helical 
flagella of E. coli are rotated at a frequency co « 100 Hz, resulting in a swimming velocity of about U « 30 /mi/s for 
a bacterium of size L « 10 fim. As a consequence, E. coli swims approximately at the speed of U « 1/30 body length 
per beat. As a difference, the flagellum of sea-urchin spermatozoon is actively oscillating at a frequency uo « 40 Hz, 
resulting in a velocity U ~ 200 /im/s for an organism of size L « 50 /am [5]. The swimming speed in this case is 
therefore about U ~ 1/10 body length per beat. How does that compare with our optimal swimmers? To answer this, 
we look at the results of Table I for optimal swimmer II. Since we have nondimensionalized times by cj _1 , we need to 
multiply the results by 2tt to obtain velocities per unit frequency. In the case e ~ 1, we find that the optimal swimmer 
with spherical body swims at about U « 1/90 body length per beat whereas the optimal swimmer with slender body 
has a velocity U ~ 1/50 body length per beat. The optimal elastic swimmers have therefore swimming performances 
which are comparable (although smaller by less than order of magnitude) with that of real microorganisms. The most 
efficient elastic swimmers (typically 0.1 — 0.4% efficiency) are also outperformed by typical swimming microorganisms 
(usually 1 — 2% of swimming efficiency). Both performances result from the exponentially decaying filament amplitude 
in the case of elastic swimmers, as compared with the large amplitude motion of real flagella. 

To conclude, we note that the results of our study could be improved upon in different ways. First, the treatment 
of the hydrodynamics of the filament using resistive-force theory is approximate and could be improved upon by using 
slender-body theory - most likely numerically. Hydrodynamic interactions between the filament and the swimmer 
body should also be included. Further improvement could be obtained by computing the swimming trajectories in 
the case of large-amplitude actuation and by including, in the case of three-dimensional actuation, the filament twist 
strains. Finally, the issue of thermal fluctuations (for the filament shape) and diffusion (for the swimmer position 
and orientation) should also be considered. Work in these directions is currently underway and will be reported in a 
future paper. 
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APPENDIX A: RELATIONSHIP BETWEEN THE BODY-FIXED AND THE LABORATORY-FIXED 

FRAMES OF REFERENCE 



In this appendix we derive the relationship between the swimming kinematics in the frame translating and rotating 
with the swimmer body and those in the laboratory frame of reference. The calculations are straightforward but have 
to be done correctly, and therefore are worth deriving. As in the main text of this paper, we denote by {e x ,e y ,e 2 } 
the cartesian coordinate system moving with the swimmer body and by {ei,e 2 ,e 3 } the cartesian coordinate system 
in the laboratory frame of reference, with ei being the average swimming direction. In the limit of e going to zero, 
we expect e x to be almost equal to ei, and both e y and e z to be almost given by a solid body rotation around ei at 
a constant rate, equal to Q x . 

Let us define the matrix 

E(t,r) = [e x e y e,] T , (Al) 
then the equation for the kinematics of the moving frame is given by 

±E = M(t)-E, (A2) 



with 



M(t) = 



( o n z (t) -n y (t) 

-n z (t) o n x | • (A3) 
\ Q y (t) -n x o 



The subtle issue which arises here is the appearance of a new (long) time scale ~ l/ft x ~ 1/e 2 , time scale over 
which the y and z axis are expected to rotate by an angle ir/2. Since the other time scale is ~ 1, then in order to 
obtain a solution for the dynamics of the body-attached frame for all times (both short and long time scales) in the 
limit of small e, we have to use the method of multiple scales. 

We then define a new time r = e 2 t, so that d/dt is formally replaced by d/dt + e 2 d/dr. In that case, Eq. (A2) 
becomes, in the multiple scale setting, 

^ + e 2 ^jE=(eM(t)+e 2 N)-E, (A4) 

where we have used the fact that both Q y and Q z are function of the short time scale and of order e, while Q x is 
constant and or order e 2 . We then look for a solution as a perturbation expansion E(t,r) = Eq + eE\ + e 2 E2 + 
and obtain 

= 0, (A5a) 

= M(t)>E , (A5b) 

= N>E + M(t)>E u (A5c) 

= N-E ± + M(t) -E 2 . (A5d) 





dE 




dt 




dE 1 




dt 


dE 2 


dE 


at 


dr 


dE 3 


dEx 


dt 


V dr 



We then can solve the system given by Eq. (A5) order by order, using the usual multiple-scales trick that terms which 
would violate the perturbation expansion hypothesis have to be set to zero, and we obtain, written in the original 
time variable for the first two terms as 

e x = ei + 3? j^e"* £ j [cos(^t)e 2 + sin(^t)e 3 ] - 3? jiO^e - ^ j [- sm(Q x t)e 2 + cos(Q x t)e 3 }, (A6a) 

ei + cos(tt x t)e2 + sin(^t)e3, (A6b) 

&z = 3? jffi^e - ^! ei — sm(Q x t)e2 + cos(^ x t)e3, (A6c) 
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where the notations introduced in the main part of this paper for the Fourier transforms have been used. We can 

now write down the relationship between velocities in the laboratory frame and velocities in the frame moving with 
the swimmer: 

U = U x e x + U y e y + U z e z , (A7) 
which means, given Eq. (A6), that we have, at leading order for each component 

U x = U x -U y ^^ih z e- it ^U z ^^ih y e- it Y (A8a) 

U 2 = UyCos(Q x t)-U z sm(n x t), (A8b) 

Us = U y sm(Q x t)-\-U z cos(Q x i). (A8c) 

Finally, denoting the averages on the short time scale by (...), we get 

(Ui) = (U x ) + ±S{U z £ll + Ufi z }, (A9a) 

U[ = U' x + )^[e- 2i \U y h z -U z h y )y (A9b) 
and since both U y and U z average to zero, we obtain (JJ 2 ) = (Us) = 0. 
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